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A LIFTED SQUARE FORMULATION FOR 
CERTIFIABLE SCHUBERT CALCULUS 


NICKOLAS HEIN, FRANK SOTTILE 

Abstract. Formulating a Schubert problem as the solutions to a system of equations 
in either Pliicker space or in the local coordinates of a Schubert cell usually involves 
more equations than variables. Using reduction to the diagonal, we previously gave a 
primal-dual formulation for Schubert problems that involved the same number of vari¬ 
ables as equations (a square formulation). Here, we give a different square formulation by 
lifting incidence conditions which typically involves fewer equations and variables. Our 
motivation is certification of numerical computation using Smale’s a-theory. 


A m X n matrix M with m > n is rank-deficient if and only if all of its n x n minors 
vanish. This occnrs if and only if there is a nonzero vector n G C"" with Mv = 0. 
There are (™) minors and each is a polynomial of degree n in the mn entries of M. In 
local coordinates for v, the second formnlation gives m bilinear eqnations in mn+n—1 
variables, and the map (M, n) i—>■ M is a bijection over an open dense set of matrices of 
rank n—1. The set of rank-deficient matrices has dimension (m-|-l)(n—1), which shows 
that the second formnlation is a complete intersection, while the first is not if m > n. The 
principle at work here is that adding extra information may simplify the description of a 
degeneracy locns. 

Schnbert varieties in the flag manifold are nniversal degeneracy loci [B] . We explain how 
to add information to a Schnbert variety to simplify its description in local coordinates. 
This formulates membership in a Schubert variety as a complete intersection of bilinear 
equations and formulates any Schubert problem as a square system of bilinear equations. 
This lifted formulation is both different from and typically significantly more efficient than 
the primal-dual square formulation of [8] , as we demonstrate in Section [3l 

Our motivation comes from numerical algebraic geometry [TB], which uses numerical 
analysis to represent and manipulate algebraic varieties on a computer. It does this by 
solving systems of polynomial equations and following solutions along curves. For nu¬ 
merical stability, low degree polynomials are preferable to high degree polynomials. More 
essential is that Smale’s a-theory [T3] enables the certification of computed solutions to 
square systems of polynomial equations [TOj, and therefore efficient square formulations of 
systems of polynomial equations are desirable. Furthermore, the estimates used in imple¬ 
mentations of a-theory simplify for bilinear systems, as explained in [HI Rem. 2.11]. Inter¬ 
estingly, formulations as square systems of bilinear equations may also aid Grobner basis 
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computations. Faugere, et al. [5] gave improved complexity bounds for zero-dimensional 
bilinear systems. 

The Schubert calculus is a well-understood, rich family of enumerative problems which 
has served as a laboratory to study new phenomena in enumerative geometry [13] . Prob¬ 
lems in Schubert calculus lead to highly-structured systems of polynomials that are chal¬ 
lenging to study. Traditional formulations of most problems in Schubert calculus are not 
complete intersections, and those which are complete intersections have far fewer solutions 
than predicted by the BKK bound |1]— this is demonstrated in Table 2 of [TB] . 

Square formulations of Schubert problems also enable the certified computation of mon- 
odromy, using either the algorithm of Beltran and Leykin [2113] or the Newton homotopies 
of Hauenstein and Liddell [9]. This will in turn enable the certified computation of Ga¬ 
lois groups [nils]. Because general degeneracy loci are pullbacks of Schubert varieties, 
these square formulations may lead to formulations of more general problems involving 
degeneracy loci as square systems of polynomials. 

In Section [T] we explain the traditional formulation of Schubert problems using Stiefel 
coordinates and determinantal equations expressing rank conditions. In Section [2] we give 
our new lifted square formulation for Schubert varieties and Schubert problems, illustrat¬ 
ing with some examples. In Section [3] we compare the efficiency of the lifted formulation 
with the primal-dual formulation of [ 8 ], demonstrating that the lifted formulation typically 
involves fewer equations and variables, and through three examples that computations us¬ 
ing it consume fewer resouorces. 

1. Determinantal formulation of Schubert problems 

The Schubert calculus involves all questions of determining the (flags of) linear sub¬ 
spaces of a vector space that have specified positions with respect to other (fixed, but 
general) flags of linear subspaces. We briefly sketch the Schubert calculus, Stiefel co¬ 
ordinates for Schubert varieties, and traditional determinantal formulations of Schubert 
problems. We work over the complex numbers for convenience and motivation from nu¬ 
merical algebraic geometry. Our formulations and main result are valid over any held, 
if we replace claims of transversality by properness (expected dimension) when the held 
has positive characteristic. This is because Kleiman’s result showing transversality of the 
intersection of general translates becomes properness in positive characteristic m- For 
the Grassmannian, we retain transversality as Vakil na proved that general translates of 
Schubert varieties in a Grassmannian intersect transversally in any characteristic. 

1.1. Schubert problems. Fix an integer n and a sequence a,: ai < ••• < < n of 

positive integers. A flag of type a, is a sequence of linear subspaces 

E, : Ea, G Ea^ G ■■■ G Ea, C C”, 

where dim Eaj = dj- A hag is complete if a, = {1, 2,..., n—1}. Given a hag E, of type 
a,, there is a list (ei,..., e^J of independent vectors such that Ea- is the linear span of 
{ei,..., 6 ( 1 ^.} for each 1 < j < s. In this case, write E, = ((ci,..., ea^))a,- The set of all 
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flags of type a, is an algebraic manifold F£(a,; n) of dimension 

S S 

(1.1) dim(a,) := — aj)(aj — aj_i) = n ■ Qg — aj{aj — aj-i ), 

j=l i=l 

where Oq := 0. When s = 1, the flag manifold ¥i{a,]n) is the Grassmannian of Oi-planes 
in C"’, Gr(ai;n), which has dimension ai(n—oi). 

The position of a flag E, of type a, with respect to a complete flag F, is the nx s array 
of nonnegative integers dim(Fj fl Ea-) for i = 1,..., n and j = 1,..., s. These positions 
are encoded by permntations w E Sn with descents in a,. For snch a permntation w, 
w{i) > tc(i+l) implies that i = a^, for some j. Write VF“* for this set of permntations. 
Given w E hF“* and a complete flag F,, we have the Schubert cell, 

(1.2) X°F. := {E, E ¥i{a.) \ dim{Fi G = #{fc < a,- | w{k) < z}} . 

The Schubert variety X^F, is the closnre of X^F, and is obtained by replacing the dimen¬ 
sion eqnality in fll.2p with an ineqnality >. This has dimension i{w) := ^{k < j \ w{k) > 
w{j)}, the nnmber of inversions of w, and thus codimension |tc| := dim(a,) — i{w). 

A Schubert problem is a list w := {wi ,..., Wr) of elements Wi E W°“ for z = 1,..., r 
satisfying Iznil |znr| = dim(a,). Given a Schubert problem w, Kleiman showed pT] 

there is an open dense subset of the product of flag manifolds consisting of r-tuples of 
flags {Fl,... ,Fl) such that the intersection 

x^,Fl n x^,Fl n n x^f: 

is transverse. Kleiman’s Theorem implies that the points of intersection lie in the corre¬ 
sponding Schubert cells—we lose nothing (for general flags) if we restrict to Schubert cells, 
and the same reasoning allows us to restrict to any dense open subset of the Schubert 
varieties. The number of points in the intersection is independent of the choice of general 
flags and this number may be determined by algorithms in the Schubert calculus. 

1.2. Determinantal formulation of a Schubert variety. Suppose that A is a set of 

nxug matrices x whose column vectors ei{x),..., ea„(x) are independent. The association 

X 3 X I—)■ ((ei(x),... ,ea,(x)}}a. =: F,{x) 

dehnes a map X -3 ¥£{a,;n). We call X Stiefel coordinates for the closure of the image 
of this map. We have Ea.{x) := span{ei(x),..., eaj{x)}, and we also write Ea.{x) for the 
n X Qj matrix whose columns are ei(x),... ,eaj{x). Whether we intend the subspace or 
the matrix will be clear from context. 

Suppose that a set A of rz x Ug matrices forms Stiefel coordinates for some subset 
X C F.^(a,; n). Let F, be a flag with a basis /i,..., /„ that forms the columns of a n x rz 
matrix. Write F^ both for the /c-dimensional subspace of the flag F, and for the n x k 
matrix with columns fi,..., fk- 

For w E hF“*, set rij{w) := < aj \ w{k) < z}. Then E, E Xy^F, if and only if 

dimFiflFa. > rij{w) 


i = l,...,n i = l,...,s. 
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Then the condition on x E X that E,{x) G X^F, is 

rank(Fi I < i + aj - rij{w) i = j = l,...,s. 

This is given by the vanishing of minors of (Fj | Eaj) of size i + aj — rij{w) + 1. These 
are polynomials in the entries of x E X. Not all such minors are needed. Even if 
redundant minors are eliminated, the number that remains will in general exceed |w|. 
This is discussed for Grassmannians in Section 1.3 of [S], where it is shown that after 
removing redundancy, |r(;| = 0 , 1 or ai = 1 , n—1 are the only cases for which this number 
of minors equals |tc|. In Section [XT] we present a typical example minimally requiring 17 
minors, but where |tc| =4. 

1.3. Stiefel coordinates for Schubert varieties. A Schubert cell X^F, has a descrip¬ 
tion in terms of bases. For a flag F,, set Fq := {0}. 

Lemma 1.1. Let w E IT“* and F, be a complete flag. Then a flag E, of type a, lies 
in X^F, if and only if there exist vectors (ei,...,eaj with eu E Fw{k) \ Fu)(fc)-i for 
k = l,...,as and E, = ((d,.. .,eafl)a.- 

Proof. Let ei,..., be such a collection of vectors with Ck E Fw{k) \ Fw{k)-i- As tc is a 
permutation, these vectors are linearly independent. If E, = {{ei,... ,eafl)a., then 

Fi n Eaj = {ck \ k < ttj and w{k) < i ), 

and so E, E X^F,. Conversely, if E, E Xf^F,, observe that if aj^i < k < aj, then the 
condition that E, lies in the Schubert cell and w E implies that 

dimF^(fc)nF„, = 1 dimF^(fc)_i n . 

For each such j and k, let be a nonzero vector that, together with F.uj{k)-i G Ea^, spans 
Fw{k) G Fdy Then E \ and F* ..., ^as))a«* ^ 

Lemma [LT] leads to the usual Stiefel coordinates for Schubert cells [H Ch. 10]. Given 
w E W°“, let AT, be the collection of n x a* matrices (xij) such that 

^w(k),k 1 foi' k 1 , . . . , ttg 

Xij = 0 if i > w{j) or i = w{k) for some k < j , 
and Xij is otherwise unconstrained. For example, here are typical matrices in for 


w = 5724613 with a. 

= (2, 5) and w = 

3652471 with 

a, = 

(2, 3, 5, 6), 




/ 

3 ^ 1,2 

Xl ,3 

3 ^ 1,4 

^ 1,5 \ 


/ 2 : 1,1 

2:1,2 

2:1,3 

2:1,4 

2:1,5 

2:1,6 \ 


X 2 ,l 

3 ^ 2,2 

1 

0 

0 


2:2,1 

2:2,2 

2:2,3 

1 

0 

0 


X 3 ,l 

X 3,2 

0 

3 ^ 3,4 

2:3,5 


1 

0 

0 

0 

0 

0 



X 4,2 

0 

1 

0 


0 

2:4,2 

2:4,3 

0 

1 

0 


1 

0 

0 

0 

0 


0 

2:5,2 

1 

0 

0 

0 


0 

X6,2 

0 

0 

1 


0 

1 

0 

0 

0 

0 


V 0 

1 

0 

0 

0 / 


V 0 

0 

0 

0 

0 

1 / 


It is an exercise to show that if x G A^,, then F,(x) G X°F,, where F, is the standard co¬ 
ordinate flag in C". Suppose that F, G X^F, and F, = ((ei,..., eaf))a. as in Lemma (HH 
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Let y he a n X as matrix with column vectors ei,..., If we reduce each column of y 
modulo those to its left, we obtain a matrix in fhu,. We summarize this discussion. 

Lemma 1.2. For any w G the set gives Stiefel coordinates for the Schubert 

variety X^F, where F, is the standard coordinate flag. The map X^^ —)■ XyjF, defined by 
X I—)■ Efix) is a bijection between X^ and the Schubert cell XfF,. 

An entry (i, j) is unconstrained for matrices in X^ when i < w{j) and there is no fc < j 
with i = w{k). As w is a permutation, there is some k > j with i = w{k). Thus the 
unconstrained entries in Xy^ correspond to inversions in the permutation w, and so we 
conclude that dimA^^, = i{w), the number of inversions in w. 

1.4. Determinantal formulation of a Schubert problem. Let w = (tci,... ,Wr) be 
a Schubert problem and suppose that Ff,..., Ff are general complete flags. Choosing a 
basis for C”, if necessary, we may assume that Ff is the standard coordinate flag. Let 
Ff,..., Ff he n X n matrices corresponding to the flags of the same name. Then, in the 
local Steifel coordinates for Xyj^Ff, the instance 

Xy^.Ff n Xy,,F^ n ■ ■ ■ n 

of the Schubert problem is given by the rank conditions 

rank ( F^ \ Ea.{x)) < i + aj - rij{wk) 

for i = 1,..., n, j = 1,..., s, and k = 2,..., r. These rank conditions are equivalent 
to the vanishing of minors of appropriate sizes of these matrices. As we discussed, this 
typically involves more equations than variables. Call this the determinantal formulation 
of the Schubert problem. 

2. Lifted square formulations for Schubert problems 

We give a new formulation for Schubert varieties as complete intersections in that the 
number of variables is equal to the sum of the dimension of the Schubert variety and 
the number of equations. These equations are bilinear when we use Stiefel coordinates 
for the flag manifold. This leads to a square formulation of any Schubert problem. In 
Subsection 12.21 we explain an improvement to this formulation. 

Fix a sequence a,: ai <■■■< as < n and let E, be a flag of type a, in C". A 
complete flag E, in C” induces complete flags on each quotient vector space EajEa^_.^ for 
j = 1,..., s. The subspaces in the induced flag on EaJEa^_.^ are 

(2.1) {fEa^r\Ek) + Ea._^)/Ea._^ foT k = 1, . . . , U . 

If w is the unique permutation in hF“* such that E, G X°F,, so that E, and E, have 
relative position w, then the subspaces fl2.ip in the flag on EajEa^_^ induced by E, are 

((-^aj Cl Eyj(j^'fi T ^for Uj —1 < fc F aj . 

(Recall that aj_i < i < k < aj implies that w{i) < w{k) and thus F^p) C F^(fc).) When 
E, = ((ei,..., ea^))a. for independent vectors ci,..., we have another complete flag 
in each quotient space EaJEaj_.^ for j = 1,..., s, whose subspaces are 

(2.2) ((cfc, Cfc+i,..., e^.) + Ea._^)/Ea._^ for a^.i < k < aj . 
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We say that (ei,..., CqJ and F, are in a,-general position if for each j = 1 ,... ,s, the two 
flags fl 2 .ip and fl 2 . 2 l) on EajEaj_i are in linear general position. That is, an intersection 
G n hr of snbspaces, one from each flag, has the expected dimension dim G + dim H — 
dim(h;„./h;„._J. 

The set of those (ei, ..., CqJ with E, = {{ei,... ,ea^)) a. that are in a,-general position 
with F, forms an open and dense snbset of those (ci,..., e^J with E, = {{ei, ... ,ea^)) a.- 
Indeed, there is a dense open snbset of the general linear gronp giving linear combinations 
of the snblist ..., Ca^ which indnce a flag on EajEa._.^ in linear general position 

with the flag indnced by F,. 


2 . 1 . Lifted square formulation. The lifted square formulation relies upon the following 
lemma. For a number k < as, define := min{aj | k < aj}, which is the smallest 

number in a, that is at least as large as k. 

Lemma 2.1. Suppose that E, = {{ci, ... ,ea^))a. is a flag of type a,, F, is a complete flag 
in a,-general position with (ci,..., e^J, and w G . Then E, G W°F, if and only if 
for each /c = 1,..., there are numbers for i < f/cla. with w{k) < w{i) such that 

( 2 . 3 ) gk . Cfc T ^ ^ , 

r^i 

w{k)<w(i) 

where gu G Fy^i^k) \ Fw{k)-i- Furthermore these numbers ak,i are the unique numbers with 
this property. 

We illustrate this lemma with two examples. 

Example 2.2. Suppose that := (61,62,63) lies in the Schubert cell 

Grassmannian Gr( 3 ; 8) and (61,62,63) is in general position with F,. Then there are 

constants 01,2, 01,3, and 02,3 such that if 

gi ■= 61 + 01^262 + 611,363, 

( 2 . 4 ) g2 := 62 + 02,363, and 

5'3 := 63, 


then gi G F3, g2 G F5, and gs e Eg = C®. 

View these now as variables and equations for membership in V35812467The linear 
forms defining the subspaces in F, give 5 + 3 + 0 = 8 equations on the vectors 61, 62, 63 and 
variables 01,2, 01,3, 02,3. As these linear forms are are general, they define a subset of codi¬ 
mension eight which when projected to the Grassmannian gives a subset of codimension 
five, which is the codimension of W35812467A,. o 

Example 2 . 3 . Suppose that E, := ((ci, 62, 63, 64))2<4 lies in the Schubert cell i^5g4Yi2368^» 
of the flag manifold F£( 2 , 4 ; 9 ) and (ci, 62, 63, 64) is in general position with F,. Then there 
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are constants ai,2, 034, 03,2, 03,4, and 04,2 such that if 


91 '■= Cl + ai^ 2&2 5 

92 '■= 62 , 

fi'a := <T 3 ,i^i + 03,262 + 63 + 03,464 5 and 

fi'a ■= 04,262 + 64, 

then gi G F5, 92 ^ Fg = C®, G F4, and g^i G F-j. As a formulation for membership in 
A59 4712368-P"#, the linear forms dehning the Fj give 4 + 0 + 5 + 2 = 11 equations on the 
vectors 61, 62, 63, 64 and hve variables o^,*. As these forms are general, they dehne a subset 
of codimension eleven that when projected to Ff'( 2 , 4 ; 9 ) gives a subset of codimension six, 
which is the codimension of ^594712368Since the membership equations {g^ G F^p), 
etc.) are linear in the variables o^,*, the hbers over points of X594712368are affine spaces. 
The equality of dimensions and surjectivity implies that the hber over a general point is 
a singleton, which is the unicity assertion in Lemma 12.11 o 


Proof of Lemma lKJi Suppose hrst that gk G F^(fc) \ Fy^(^k)-i where ^fi,..., ga^ are dehned 
using fl 2 . 3 p for some constants ak,i. Then F, = ((^fi,... ,ga^))a., as the expressions ( 12 . 31 ) 
are unitriangular. Lemma [ 1.11 then implies that F, G XfF,. 

For the other direction, we use induction on j to construct unique constants such 
that the vector gk dehned by ( 12 . 31 ) satishes gu G F^(fc) \ F^(^k)-i for k < aj. We will 
suppose that that for each k < aj_i there are unique constants for i < [fc]a. with 
w{k) < w{i) such that if gk is the linear combination ( 12 . 3 p . then gk G F^(fc) \ F^(fc)_i, and 
use this to obtain the constants 0^,^ for aj_i < k < aj. This is no assumption in the base 
case (j = 1) of this construction. 

By our assumption on (ei,..., CaJ and F,, the two hags in Fa^./Fa^_^, 


Faj_i ^ + (fiaf) ^ ^ + (6(j^._j^-|_ 2, ..., 60^) ^ ’ aud 

Faj-i $ A {Fw{aj-i+l) n Eaj) C • • • C Eaj_-^ + (F^(a^._i) H Eaj) F F, 
are opposite. In particular, for any aj-i < k,i < aj, we have that 

(2.5) + iF^(^k) (d Fa^)^ n (^Eaj_^ + (6j, . . . , 6(J^.)^ 


has dimension max( 0 ,/c+ 1 —z) modulo This implies that there are constants ak,e 

ioT k < £ < Oj and an element e G Fa^_^ such that the sum 


( 2 . 6 ) 


6 fc + Oik/Ci + 

e=k+i 


e 


lies in F^(fc). In fact, the sum (12.6p lies in F^(fc) \ F^(fc)_i. Indeed, as F, G X°F,, we have 
that F^(fc)_i n Fa. C F^(fc) n Fa., and so the dimension of (12.5p drops if we replace F^(fc) 
by Fw{k)-i- This also implies that the numbers ak,e are unique. 

The element e G Eaj_i is some linear combination of 6i,...,eaj_j and thus also of 
9 i) ■ ■ ■) 9aj.i- Since gi G F^p), those gi with w{i) < w{k) are not needed for the sum (12.61) 
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to lie in and thus there are constants ( 3 i for i < aj-i with w{k) < w{i) such that 

aj 

( 2 . 7 ) Qk := Cfc + ^ ak,iee + ^ ( 5 igi 

i=k-\-l i<aj 

w{k)<w{i) 

lies in \ Fy^(^k)-i- As each gi in the second sum lies in \ the constants 

Pi are unique. To obtain the expression fl 2 . 3 p for gk hrst use the formula fl 2 . 3 p for each gi 
appearing in fl 2 . 7 p to rewrite the second sum as a linear combination of for ^ < aj_i 
with w{k) < w{i) < w{i), and then use that w E hh“* to see that {fc+ 1 ,..., ajj is the set 
of z in the interval (aj-i,aj] with w(k) < w(i). The unicity of the constants ak^i follows 
from that of the constants ak^i and Pi, and our induction hypothesis. □ 

Remark 2 . 4 . Lemma o leads to a square formulation for membership in X^F, for flags 
in F£(a,;n) as follows. 

( 1 ) Pick Stiefel coordinates Xa, for ¥i{a,;n). For x E Xa,, we have the partial flag, 
F,(x) = ((ei(x ),.. .,ea,(x)}}a. E F£(a,;n). 

(2) Choose lifting coordinates 

(2.8) a = {ak,i j k = 1 ,..., ag, i < \k~\a. with w{i) > w{k)} 
and form the vectors 

gk{x,a) := ek{x) + ^ ak^iC^x), 

w(k)<w{i) 

for /c = 1 ,..., as- 

( 3 ) Given independent linear forms /i,..., /„ such that Fj is dehned by the vanishing 
of /j+i ,..., fn, our equations for E,{x) E X^F, are 

fjigpx, a)) = 0 for z = 1 ,..., and j > w{i). 

These equations are bilinear in the sets of variables x E Xa, and a. o 

Definition 2 . 5 . Call the formulation for membership in X^F, for flags in F£(a,;n) of 
Remark 12.41 the lifted formulation for a Schubert variety. Write a{w) for the set of lifting 
coordinates fl2.8ll and |a(zu)| for the number of these coordinates, which is 

ds 

( 2 . 9 ) |a(zu)| = ^#{z < I ^(z) > zu(fc)}. 

k=l 

Theorem 2 . 6 . The lifted formulation for membership in X^jF, C F.^(a,) is a complete 
intersection. 


Proof. We must show that dimX^F, equals the number of variables minus the number of 
equations. The number of equations is the sum of codimensions of the Fw{k) for k < Og, 


ds 

''^^n — w{k) = 

k=l 


k=l 


( 2 . 10 ) 


n ■ a 












A LIFTED SQUARE EORMULATION FOR CERTIFIABLE SCHUBERT CALCULUS 


9 


The number of variables is the dimension of Ff'(a,;n), as calculated in fll.ip 

S 

( 2 . 11 ) dim(a,) = n • — aj-i)aj , 

i=i 

where Oq = 0 , plus the number |a(tc)| of the variables We rewrite fl 2 . 9 p as 

CLs 

k=l 

CLs 

k=l 

The hrst equality uses that if aj < i < k < Oj+i, then w{i) < w{k), as w E 1 T“*. We 
rewrite this as 

S 

(2.12) ( — aj_i)aj^ — < k < ag \ w{i) < w{k)} . 

i=i 

Using that w E W°“, the linear combination fl 2 . 12 |) + f| 2 . 1 ip — fl 2 . 10 p becomes 

ds 

< k < ttg \ w{i) < w{k)} 

k=l 

= #{i < k I w{i) > w{k)} = i{w) = dimX^F,, 

which completes the proof. □ 

Remark 2 . 7 . The lifted formulation of Remark 12.41 for a Schubert variety leads to a 
square formulation for Schubert problems, following Subsection 11.41 Suppose that w : = 
(tci,..., tCr) is a Schubert problem on F£(a,;n). Let be general flags and 

consider the intersection of Schubert varieties 

( 2 . 13 ) X^,F^ n X^,F^ n ■ ■ • n 

Assume that F^ is the standard coordinate flag and use Steifel coordinates for the 
Schubert cell X^^F^ to formulate the intersection fl 2 . 13 l) . Replacing the determinantal 
rank conditions for membership in each Schubert variety X^jj^F^, ■ ■ ■, X^^Fl by the lifted 
square formulation gives the lifted formulation for the Schubert problem w. It uses 

£(wi) + |a(tC 2 )| H-h |a(tCr)| 

variables and bilinear equations. o 

Since the intersection fl 2 . 13 p is transverse, it is zero-dimensional (or empty). This gives 
the following corollary to Theorem 12.61 

Corollary 2 . 8 . The lifted formulation for membership in the intersection (I 2 . 13 p is a 
complete intersection in the local coordinates . 
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Remark 2.9. For Grassmannians, there are Stiefel coordinates parametrizing the intersec¬ 
tion n [HI § 3.1]. These involve dim(Gr(ai;n)) — |tci| — \w2\ = — 1102! 

variables and lead to a lifted formulation of fl2.13p using 

t{wi) - |tC 2 | + |a('ii'3)| H-1- \<y{Wr)\ 

variables and bilinear equations. This presents fl2.13p as a complete intersection using 
\w 2 \ + \a{w 2 )\ fewer equations and variables than the formulation of Gorollary 12.81 o 

2.2. Reduced lifted formulation. We introduce an improvement to the lifted square 
formulation, motivating it through three examples. 

Example 2.10. Gonsider the lifted formulation for the Schubert variety in Gr(3; 8 ) 

where w = 45812367. Suppose that = ( 61 , 62 , 63 ) where ( 61 , 62 , 63 ) come from Steifel 
coordinates X for Gr(3; 8 ) and involve 15 variables. The lifted formulation uses three 
new variables 01 , 2 , «i, 3 , and 02,3 as in Example 12.21 and we form the vectors gi,g 2 ^gz as 
in fl2.4l) . Then E, G X^F, if and only if gi E F 4 , g 2 E F 5 , and g^ E Fg, giving seven 
equations in 15 -|- 3 variables to dehne the codimension four Schubert variety X^F,. 

It suffices to only require that gi and g 2 lie in F 5 , for then some linear combination of 
the two will he in F 4 . This dispenses with one equation. Having done this, we may also 
dispense with the variable ai^ 2 , and thereby obtain a reduction of one variable and one 
equation. Specihcally, suppose that 

gi '■= Cl -7 ai,363 , 

g 2 := 62 + 0 ^ 2,363 , and 

gs ■= 63 . 

Then F 3 G X^F, if gi,g 2 E F 5 . This gives six equations in 15-1-2 variables to define X^F, 
in the Steifel coordinates X for Gr(3;8). o 

Example 2.11. A similar reduction is possible in Example 12.31 The requirement that 
(73 G F 4 may be relaxed to g^ G F 5 , for then some linear combination of g^ and gi lies in 
F 4 . This removes one bilinear equation, and we may dispense with 613 4 . o 

Example 2.12. Now consider the lifted formulation for Xy^F, with w = 358 47126, which 
has codimension six in the 21-dimensional flag manifold F£(3, 5; 8 ). Suppose that 61 ,..., 65 
are the column vectors from the Steifel coordinates <T' 3<5 for F£(3, 5; 8 ) and set E, := 
(( 61 ,..., 65 )) 3 < 5 . The lifted formulation uses seven variables ak,i to form the vectors 


9 i = 

61 -|- ttl,262 + 6^1^363 , 


g 2 = 

62 + 0^2,363 5 


93 = 

es, 


94 = 

0^4,262 + 6^4^363 -1- 64 -f- 6^4^565 

9 b = 

as, 363 + 

es 


which are required to he in the subspaces 

S'! ^ -^3 ) 5'2 G F5 , (73 G Fg , g4 E F4 , and g^ E F-j. 
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This gives 5 + 3 + 0 + 4 + 1 = 13 bilinear equations in |a,| + |a(ta)| = 21 + 7 variables to 
define XyjF,. 

Some of these conditions and variables are redundant. All that is needed is that Qi & F^, 
92,94 £ As and 93,95 G Ag. The first three give 5 + 3 + 3 = 11 bilinear equations and the 
last two give none. Similarly, the variables 04 ^ 2 , and 05^3 are not needed. Thus X^F, has 
a formulation involving five new variables and eleven bilinear equations. o 

The reduction in these examples was possible when for some k < a.^ there was a number 
m such that the consecutive values w{k)+l,..., w{k)+m for the permutation w occured at 
positions i < \k~\a,- If A, • • • Am < \k]a, are the positions such that w{ij) = w{k)+j, then 
the condition of Lemma 12+1 that Qk G A^(fc) may be replaced by Qk G Fyj{^k)+m = Fw{im), for 
there is some linear combination of the vectors gk, Qi^,..., that lies in A^(fc). Likewise, 
the variables ak^^, ■ ■ ■, C(k,im are not needed. 

We formalize this. Consider vectors ei(a;),..., ea^{x) coming from Steifel coordinates X 
for some subset X of F.^(a,; n). For each /c = 1,..., a^, let f3k be the set of indeterminates 

(2.14) 13k := {(3k,i I i < \k]a. and 3 j > \k]a. with w{k) < w{j) < w{i)} , 
and set 

(2.15) Qk = 9k{x,(3) := ek{x) + '^Pk,iei{x) . 

i 

Write (3 = (3{w) = Uk(3k for the set of all these indeterminates and |/3(tc)| for the number 
of indeterminates in (3. 

For a complete flag F,, the reduced lifted formulation for membership in X^F, in the 
Steifel coordinates X uses the additional variables I3{w) to form the expressions fl2.15p . 
and has the equations given by the membership requirements 

9 kix, /3) G Fyj(^k) 4 -m{k) , fur k 1, . . . , Clg , 

where m{k) is the largest number m such that the consecutive values w{k)+l ,..., w{k)+m 
for the permutation w occur at positions i < [/c]a. • 

The results in Subsection 12.11 hold mutatis mutandis for this reduced lifted formulation 
of Schubert varieties and Schubert problems and are omitted. 

3. Comparison with the primal-dual square formulation 

We compare the efficiency of this lifted formulation to the primal-dual formulation 
of [H] . Both involve added variables and bilinear equations in local Steifel coordinates. We 
first compare these formulations to the determinatal formulation of a particular Schubert 
variety. Next, we determine which of the two uses fewer added variables for each Schubert 
variety on a flag manifold in C®, and then compare their computational efficiency for 
solving three Schubert problems, including two from [ 8 ]. We almost always observe a gain 
in efficiency for the lifted formulation over the primal-dual formulation. 

We may take advantage of whichever formulation is most efficient for a given Schubert 
variety, for they are compatible. That is, one may construct a hybrid system of equations 
for the intersection fl2.13p using a lifted formulation to determine membership in some 
of the Schubert varieties and a primal-dual formulation to determine membership in the 
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others. Whenever |ty| = 1, the determinantal formulation for membership in the hyper¬ 
surface Schubert variety is a single determinant in Steifel coordinates, so there is 

no need to use an alternative formulation to obtain a square system. In what follows, we 
will always use the determinantal formulation when |tc| = 1 . 

3.1. Comparison of three formulations. We compare the three formulations, deter¬ 
minantal, primal-dual, and lifted, for membership in the Schubert variety X 34781256 X, in 
Gr(4, 8 ). Let X be Steifel coordinates for Gr(4, 8 ), which is a set of 8 x 4 matrices of rank 
4 and F, be a flag in C®. The Schubert variety X 34781256 F, consists of those 4-planes H 
that meet the hxed 4-plane F 4 in a subspace of dimension 2. 

If F 4 is represented as the column space of a 8 x 4 matrix, then a 4-plane H from X 
lies in X 34781256 X, if and only if 

rank {H F4) < 6 . 

A priori^ each of the 7x7 minors of this matrix must vanish for a total of 64 quartic 
equations in the entries of 77 G X. In [51 § 1.3] the Plucker embedding of the Grassmannian 
is used to give a smaller set of equations which are linear combinations of the maximal 
4x4 minors of X. The dimension of the linear span of such equations is the cardinality 
of the set {p G | p ^ 3478} of increasing sequences p of length 4 from [ 8 ] = {1,..., 8 } 
where one of the inequalities pi < 3, p 2 < 4 ,p 3 < 7, or p 4 < 8 does not hold. There are 
seventeen such sequences 

5678,4678 , 3678,4578,2678 , 3578,4568 , 1678 , 2578 , 

3568,4567, 1578, 2568, 3567, 1568, 2567, 1567, 
so that in Steifel coordinates, X 34781256 X, is dehned by 17 equations. 

The primal-dual formulation uses a variant of the classical reduction to the diagonal. 
Gonsider the map T on G(4, 8 ) which sends a linear subspace H to its annihilator, 

This is an isomorphism in which T (X^F,) = X^xF,-*-, where F,-*“ is the flag of linear forms 
annihilating the linear subspaces in F, and w-^ = wqwwq, where wo{i) = n+l—i. 

To understand this in Steifel coordinates, pick a basis corresponding to the rows of a 
matrix whose dual basis corresponds to the columns. The dual Schubert variety X^,xF,-*- 
has Steifel coordinates X^±, where we send K G X^i^ to the row span of where 

4) is the matrix whose hrst i columns span Fj. In this formulation, the intersection of 
X 34781256 -^» with the set parametrized by X is the intersection of the graph of T with the 
product X 3478 i 256 iX,-*- X X. Since 3478 1256“‘“ = 34781256, the primal-dual formulation 
uses the coordinates X 34781256 x X with the equations 

= 04X4, 

which state that the four-plane annihilates H. This involves 12 = dim X 34781256 

new coordinates and 16 bilinear equations, which are the entries of the matrix . 

Finally, the lifted formulation uses the local coordinates 

y ^ ^ ^ 1.1 2 / 1.2 V 
V 0 1 1 / 2,1 1 / 2,2 / 

from Gr(2,4): For Y G y and 77 G X, the 8x2 matrix HY is a two-plane in 77. 
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If 01 ,..., 04 are the equations the dehne F 4 , then the lifted formulation for the inter¬ 
section of X 3478 i 256 -^» with the set parametrized by X uses the coordinates X x y and 
the equations 

(j)i{HY) =0 for i 

This involves 4 = dim y new coordinates and 8 bilinear equations (linear in the the entries 
of YH) as (j)i{HY) = 0 gives two equations, one for each column in HY. 

3.2. Added variables for Schubert varieties on flag manifolds in C®. The square 
primal-dual formulation of a Schubert variety on the flag manifold [8] uses that every 
flag E, in C"' has an annihilating dual flag in the dual space to C”. If E, has type 

a,, then E^ has type := {n—aj \ aj G a,}. This duality gives an isomorphism 
T: F£(a,;n) —)■ F£(a^;n) with Y{Xy^E,) = X^±E^ (we refer to Section 4 of [8] where 
vj-^ is dehned). A variant of the classical reduction to the diagonal allows us to formulate 
membership of a flag E, in X^E, by parametrizing X^±E^, using i{w) new variables. 

As explained in [SI Rem. 4.10], sometimes membership of a flag E, of type a, in a 
Schubert variety X^E, is equivalent to the membership of a projection vr(i?,) in the 
projected Schubert variety 7i{Xu,E,) = X^E,, where tt: F£(a,; n) —)■ ¥i{b,; n) is the natural 
projection and 6, C a,. When this occurs, the primal-dual formulation uses fewer, 
new variables. This is the reduced primal-dual formulation. 

For every Schubert variety X^E, on a flag manifold ¥£{a,] 9) with 1 < |t(;| < ^ dim(a,) 
we compared the numbers of new variables needed in the two formulations. The restriction 
1 < ltd I is because the determinantal formulation when |t(;| = 1 is already a complete 
intersection. The restriction |t(;| < |dim(a,) is because, as in Remark 12.71 if |t<;| > 
|dim(a,), then we would work in local Steifel coordinates A0 for the Schubert variety 
XjjjE, in any Schubert problem involving w (and any Schubert problem has at most one 
permutation satisfying this inequality). 

There are 3,395,742 such Schubert varieties in the 256 flag manifolds F.^(a,;9). We 
compared the reduced lifted formulation of Subsection 12.21 with the reduced primal-dual 
formulation for all these Schubert varieties. In 141,256 (4.160%) the primal-dual formu¬ 
lation used fewer new variables, in 3,161, 233 (93.094%) the lifted formulation used fewer 
new variables, and in 93, 253 (2.746%) the two were tied. 

This overstates the efficiency of the primal-dual formulation. For example, in only 7 
of 1725 relevant Schubert varieties in ¥i{2, 3, 5; 9) did the reduced primal-dual formula¬ 
tion involve fewer additional variables. In contrast, on the isomorphic dual flag variety 
F£(4,6, 7;9) in 124 out of 1725 relevant Schubert varieties the reduced primal-dual for¬ 
mulation involved fewer variables. 

To gain an idea of how this might be exploited, we determined which of each pair of dual 
flag manifolds F£(a,; 9) and F£(a0; 9) was more favorable for the reduced lifted formulation 
of its Schubert varieties. We redid our computation comparing the two formulations, but 
restricted it to those flag manifolds Ff'(a,;9) where Ff'(a,;9) was more favorable for the 
reduced lifted formulation than F£(a0;9). This is a fair restriction, for the number of 
additional variables in the reduced primal-dual formulation is the same for a Schubert 
variety and for its dual, but may be different for the reduced lifted formulations. 
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Redoing the computation, there were 1,877,752 Schubert varieties, as we only con¬ 
sidered one of each dual pair of flag manifolds. In 53, 698 (2.860%) the primal-dual 
formulation used fewer new variables, in 1,784,646 (95.04%) the lifted formulation used 
fewer new variables, and in 39,408 (2.099%) the two were tied. 

The reduced lifted formulation is always better for the Grassmannian Gr(fc, n) than for 
its dual Gr(n—fc,n) when 2k < n. 

Lemma 3.1. If 2k < n, then the reduced lifted fromulation always uses fewer vari¬ 
ables than the primal-dual formulation for Schubert varieties X^F, in the Grassmannian 
Gr(/c,n) with |tc| < ^{k{n—k). 

Proof. The original lifted formulation for Schubert varieties in the Grassmannian Gr(/c, n) 
used ( 2 ) additional variables, while the primal dual formulation for uses i{w) = 

k{n—k) — |t(;| variables. The lemma follows as £{w) > ^k{n—k) > Ik{k—1). □ 

Remark 3.2. The Grassmannian Gr(/c, n) has a more efficient primal-dual formulation that 
uses the Steifel coordinates of Remark 12.91 for the intersection of two Schubert varieties. 
This involves k{n—k) — |wi| — \w 2 \ new variables, while the lifted formulation uses k{k—l) 
new variables to formulate membership in two Schubert varieties. The lifted formulation 
is more efficient when 

k{n—k) — \wi\ — \w 2 \ > k{k — 1). 

Since we may assume that |wi| -|- \w 2 \ < \k{n—k), the lifted formulation is always more 
erhcient when k < (n-|-2)/3 for then 

k{n—k) — \wi\ — \w 2 \ > ^k{n—k) > ^k{2k — 2) = k{k—l). o 

3.3. Computational time and resources. We computed instances of three Schubert 
problems using the (reduced) lifted formulation. Two were computed using a primal-dual 
formulation in [8] , and the third is a problem with many more solutions. In all, the lifted 
formulation used fewer variables and less computational resources. 

Example 3.3. Gonsider the Schubert problem in Gr(3; 9) given by the permutations 

tai,...,tC 4 = 489123567 and tC 5 ,...,tCio = 689123457. 

This has 437 solutions and asks for the 3-planes in C® which nontrivially meet four given 
4-planes and six given 6-planes. The classical formulation of the intersection fl2.13p in 
Stiefel coordinates for X^^^Ff is a system of 12 variables, 20 independent linear 

combinations of cubic minors and six cubic determinants. 

The square primal-dual formulation with similar coordinates involves 24 variables, 18 
bilinear equations, and six cubic determinants. The determinants correspond to the con¬ 
ditions ws,..., wio as |689123457| = 1. In [8] we used Bertini [T] to solve an instance 
of this Schubert problem given by random real flags. This computation consumed 20.37 
gigaHertz-hours to calculate 437 approximate solutions. We then used rational arithmetic 
in alphaGertihed in to certify the solutions, which used 2.00 gigaHz-hours. 

We formulate this Schubert problem using the lifted formulation. We use Stiefel coor¬ 
dinates for Xuj^Ff n X^j^Ff which use dim(Gr(3;9)) — |tci| — \w 2 \ = 18 — 3 — 3 = 12 
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variables. The reduced lifted formulations of for i = 3,4 require a total of 

|/d(iT 3 )| + |/5(ta4)| =2 + 2 = 4 new variables and 2 • 5 = 10 bilinear equations. As in 
the primal-dual formulation, we formulate membership in the six remaining hypersurface 
Schubert varieties using six cubic determinants. The result is a system of 12 + 4 = 16 vari¬ 
ables and 10 + 6 = 16 equations. To compare with the primal-dual formulation, we solved 
a random instance using regeneration with the same variables, hardware, software, and 
software version. The lifted formulation of 16 variables and equations was a signihcant 
improvement, using only 4.75 gigaHertz-hours to calculate 437 approximate solutions. 
The output suggests 107 of the solutions are real, while the rest are non-real. We used 
seven processors in parallel, but many more could be efficiently used as the regeneration 
tracked up to 2, 265 paths in one step. 

Certihcation time was also signihcantly improved by using this formulation. A 33.54 
gigaHertz-minute computation in alphaCertihed uni using rational arithmetic verihed that 
the 437 points in the output are indeed approximate solutions and that the corresponding 
exact solutions are distinct. This computation also proved the reality for 107 of the exact 
solutions. o 

We compare the primal-dual and lifted formulations in a more general flag manifold. 
Example 3.4. Consider the Schubert problem with 128 solutions in Ff'(2,4, 5;8) given by 


Wi 

= 48 573126, 

W2,W3 

= 78 45 3126, 

W4,W5 

= 68 574123, 

Wg,W7,W8 

= 78 46 5123, and 

Wg 

= 4738 5126. 


Applying all improvements given in [5] produced a primal-dual formulation with 41 vari¬ 
ables, 36 bilinear equations, two quadratic determinantal equations corresponding to the 
hypersurface conditions w^, and three quartic determinantal equations from wq, tcy, ws- 
This square system corresponding to a random choice of nine real flags took 2.95 gigaHertz- 
days of processing power to solve and 1.78 gigaHertz-hours to certify. 

We analyze this Schubert problem using a reduced lifted formulation in the Stiefel 
coordinates consisting of i{wg) = 16 variables. The reduced lifted formulations of 
Xyj.Fl for i = 1,2,3 add |/3(tci)| + \(5{w2)\ + \/3{w3)\ = 5 + 6 + 6 = 17 new variables 
and 10 + 9 + 9 = 28 bilinear equations. As in the primal-dual formulation, we formulate 
XwiFl for i = 4,..., 8 using two quadrtic and three quartic determinants. The reduced 
lifting uses 16 + 17 = 33 variables and 28 + 2 + 3 = 33 equations. As in Example 13.31 we 
compare this with the primal-dual formulation using the tools which were utilized in [8]. 
To facilitate certihcation, we computed approximate solutions with two extra digits of 
precision compared to our computation in [8]. 

With these tighter parameters, we still observed an improvement in efficiency when 
solving a system with the new formulation of 33 variables and equations. This used 1.13 
gigaHertz-days of computing; less than half the power consumed by the similar instance 
using the primal-dual formulation. The output was 128 approximate solutions, of which 
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42 appeared to be real. Again, we used alphaCertified with rational arithmetic to certify 
the approximate solutions, verify they correspond to distinct solutions, and prove that 42 
exact solutions are real. Certihcation required 1.68 gigaHertz-hours of processor power, 
marginally less than certihcation for the similar instance we solved via a primal-dual 
formulation. 

The initial computation used six processors in parallel, but many more could be effi¬ 
ciently used as the regeneration tracked up to 708 paths in one step. Certihcation could 
have efficiently used 128 processors. o 

We formulated and solved a higher-degree problem in a Grassmannian. 

Example 3.5. Consider the Schubert problem with 28,490 solutions in Gr(3; 10) given by 

Wi, tC 2 , lys = 5910 1234678 and tC 4 ,..., Wis = 7910 1234568 . 

This asks for the 3-planes in that nontrivially meet three given 5-planes and twelve 
given 7-planes. In the determinantal formulation, we parametrize X^^Fl fl using 

dim(Gr(3; 10)) — |tci| — |w 2 | = 21 — 3 — 3 = 15 variables, and membership in is given 

by the vanishing of ten independent linear combinations of cubic minors. Including the 
cubic determinants for X^.Fl for i = 4,..., 15 uses 15 variables and 22 cubic equations. 

The primal-dual formulation uses 33 variables, 21 bilinear equations, and twelve cubic 
determinants. The lifted formulation begins with Steifel coordinates involving 15 variables 
that parametrize X^^Fl fl X^j^F^. The reduced lifted formulation for X^^F^ uses hve 
bilinear equations and adds |/3(r(;3)| = 2 variables for a total of 17 variables. The twelve 
hypersurface conditions W 4 ,..., W 15 are each given by a single cubic determinant for a 
total of 5 -|- 12 = 17 equations. The only diherence is for X^^F^ which use 18 variables 
and 21 bilinear equations with the primal-dual formulation but only two variables and 
hve bilinear equations for the lifted formulation. 

We chose 15 random real hags and solved the corresponding instance of the Schubert 
problem using 1.71 gigaHertz-months of processing power to apply regeneration in Bertini 
V. 1.4 and 4.00 gigaHertz-hours of power to apply four Newton iterations to the output 
using alphaCertihed. This produced 28,490 approximate solutions, and 1,436 appeared 
to be real. The main calculation in Bertini used 8 processors in parallel, but many more 
could be used efficiently as the regeneration tracked up to 148,161 paths in one step. 

Due to the size of the output, we soft certihed our results using hoating-point arithmetic 
in alphaCertihed with 192-bit precision. This heuristically verihed that the 28,490 points 
are approximate solutions, that they correspond to distinct solutions, and that 1,436 of 
them correspond to real solutions. This computation consumed 47.15 gigaHertz-minutes 
of processing power. A rigorous computation using rational arithmetic, but only for the 
1,436 apparently real solutions, used 1.80 gigaHertz-days and proved that 1,436 points 
in the output are approximate solutions corresponding to distinct real solutions. o 

We give additional details for the computations and comparisons in Examples 13.3113.41 
and 13.51 at the following site. 

http://www.unk.edu/academics/math/_files/square.html 
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